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C^ \ Abstract: Using a recently developed effective field theory for the interactions of nucleons 

at non-relativistic energies, we calculate the rate for the fusion process p + p —f d + e'^ + v^ to 

leading order in the momentum expansion. Coulomb effects are included non-perturbatively in a 

^S| I systematic way. The resulting rate is independent of specific models for the strong interactions at 

short distances and is in agreement with the standard result in the zero-range approximation. 

> 

\^ [ The first step in the different nuclear processes in the Sun which generate the observed 

^D ' luminosity is proton-proton fusion p+p ^ d + e"*" + z^e • It was explained more than sixty 

S— : ■ years ago by Bethe and Critchfield||2| when nuclear physics was still in its infancy. When 

Q^ ! the field had more matured, it was reconsidered in the light of more modern developments 

^N I by SalpeterO] and later by Bahcall and MayH. But in spite of the enormous progress in 

nuclear physics during this time, the methods and approximations made in these different 
calculations were essentially the same. The obtained accuracy in the obtained fusion rate 
Q I was just a few percent. Including higher order electromagnetic and strong corrections 

the uncertainty in the rate is now around one percentQl^. This is very impressive for a 
strongly interacting process at low energies very ordinary perturbation theory cannot be 
used. 
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H I In the light of the importance this fundamental process plays in connection with the 

solar neutrino problem and possible neutrino oscillations |l|, it is natural to reconsider this 
process from the point of view of modern quantum field theory instead of the old potential 
models used previously. A first attempt in this direction was made by Ivanov et al.||^]. 
They obtained then a result which was significantly different from the standard result based 
upon potential models. Subsequently it was pointed out by Bahcall and Kamionowski|^] 
that their effective nuclear interaction was not consistent with what is known about proton- 
proton scattering at low energies where Coulomb effects are important. 

The approach of Ivanov et al.[^ is based upon relativistic field theory and should in 
principle yield reliable results. But it is well known that in particular for bound states 
like the deuteron it is very difficult to use consistently a relativistic formulation. Also the 
uncertain nuclear physics part of the fusion process under consideration takes place at 
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low energies and should therefore be described within a non-relativistic framework. Then 
all the large-momentum degrees of freedom are integrated out and one is left with an 
effective theory involving only the physically important field variables. The underlying, 
relativistic interactions are then replaced by non-renormalizable local interactions with 
coupling constants which must be determined from experiments at low energies. In such a 
non-relativistic theory the important Coulomb effects can also systematically be included. 

A step in this direction has recently been taken by Park et al. using chiral perturbation 
theory in the low-energy limit [^ . They obtain results in very good agreement with previous 
potential calculations which is not so surprising since they make use of phenomenological 
nucleon wavefunctions which fit low-energy scattering data very well. A more fundamental 
approach to the same problem has recently been formulated by Kaplan, Savage and Wise in 
terms of an effective theory for non-relativistic nucleons|I^ . It involves a few basic coupling 
constants which have been determined from nucleon scattering data at low energies. With 
no more free parameters it can then be used to make predictions for the deuteron form 
factor and quadrupole moment [pT|] , deuteron polarizabilities and Compton scattering on 
deuteronsO]. The obtained results are in good agreement with experimental data. In 
this approach higher order corrections can also be derived in a systematic way. Going to 
higher energies, the effects of pions must be included using the established counting rules. 
These will cause the well-known D-mixture into the deuteron wavefunction. 

Most recently, proton-neutron fusion p + n — > d + 7 has been calculated in this theory 
by Savage, Scaldeferri and Wise[^] including the effects of virtual pions. When the process 
is taking place at very low energies, one can omit the effects of pions and replace them by 
slightly different couplings of the nucleons alone. From a field-theoretic point of view this 
process is very similar to p + p — > d-\- e'^ + v^. The main difference is the strong Coulomb 
effects which is present in the proton-proton channel. These have now been calculated and 
shown to give both a scattering length and an effective range for pp low-energy scattering 
in agreement with data|15]. Based upon these results, we will here derive the rate for the 



corresponding pp fusion reaction. This same process is also being considered by Savage and 
Wise[p^. At this stage we are only interested in the dominant contributions to the fusion 
rate in order to provide a rough comparison with results based upon potential models. 
This corresponds to the zero-range approximation in the potential model approach. In 
light of previous applications, we expect the resulting accuracy to be around 20% or 
better in the resulting rate. In this leading order approximation we ignore effective-range 
corrections, D-wave admixture, vacuum polarization, two-body current interactions and 
unknown counterterms. 

The strong interactions among the nucleons are now described by the effective La- 
grangian of Kaplan, Savage and Wisejl^]. Denoting the nucleon field of mass M by N[x) 
and including only the lowest order interaction term in the S'-channel, it can be written 
as 

£ = iVt U + Zi J N - CoiN^UN) ■ (iV^niV)t (1) 



where the IIj are projection operators into specific spin and isospin states. More specif- 
ically, for spin-singlet interactions IIj = (T2T2Tj/\/8 while for spin-triplet interactions 
IIj = cr20"ir2/v8. Calculating now the scattering amplitude for two nucleons, one finds that 



the coupling constant Cq is determined by the scattering length unn in this channel [ 10 1. 



For neutron-proton interactions in the spin-triplet channel there is pole in the corre- 
sponding Green's function corresponding to the deuteron. The residue of the pole gives 
the bound state wavefunction. With only the lowest order contact interaction in (|l|) it is 
found to be of the form 



in momentum space. This corresponds to the standard Yukawa form in coordinate space 
where I/7 represents the size of the deuteron. It determines the value of the corresponding 



coupling constant Cq in this channel|ll] 



In the absence of strong interactions, the incoming proton-proton state with center- 
of-mass momentum p is given by the Coulomb wavefunction 1 16] 



-, 00 
V'p(r) = - E(2^ + l)i'e''''Fe{p)P,icos{e) (3) 

Here p = pr and ai = argr(l + i + irj) is the Coulomb phaseshift where the parameter 
r] = aM/2p characterizes the strength of the Coulomb interaction. At low energies only 
the S-wave will contribute. It is given in terms of the confluent hypergeometric or Kummer 
function M(a, b; z) as 

Fo(p) = Cr,pe-'PM{1 - ir,, 2; 2ip) (4) 

where the normalization factor Crj = e~'^^'^|r(l + ir])\. The probability IV'plO)!^ to find 
the two protons at the same point is thus equal to the Sommerfeld fact or pT 



At very low energies when r] gets large it becomes exponentially small and is the dominant 
effect in the fusion reaction. 

The available energy in the process is set by the neutron-proton mass difference and 
the deuteron binding energy B = 2.225 MeV. It corresponds to a momentum 7 = VBM 
of the bound nucleons equal to 45.71 MeV. The temperature in the core of the Sun is 
approximately 15 x 10^ K which corresponds to an average proton momentum around 
p = 1.5 MeV. The kinetic energy of the lepton pair will therefore be much smaller than 7 
and it is a very good approximation to just ignore it. In the following the weak current is 
therefore assumed to carry zero momentum and the evaluation of the different Feynman 
diagrams will much simplify. 



To lowest order in the strong interaction between the protons, the transition matrix 
Tfi is given by the first diagram in Fig.l. After being hit by the weak current, the proton- 
proton system is transformed into a bound deuteron. The value of the diagram is thus 
seen to be 



A{p) = v/8^ 
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1 



(27r)3 k2 + 72 



^p(k) 



(6) 



where ipp (k) is the Fourier transformed wavefunction of the incoming proton-proton scat- 
tering state. To first order in the four-proton ^Si coupling Co we have the diagram in Fig. 
lb. Within the loop the protons move in the Coulomb field of each other. This motion is 







Figure 1: Feynman diagrams which contribute in leading order to proton-proton fusion. The solid 
lines are nucleons and the wiggly line represents the weak current. Between the protons there are 
exchanged Coulomb photons represented by the shaded blobs while the crossed circles represents 
the deuteron wavefunctions. 



described by the standard Coulomb propagator Gc{E) = 1/{E — Hq — Vq + ie) where Hq 
is the free, non-relativistic Hamiltonian and Vc is the Coulomb potential. In momentum 
space it takes the form 



GciE;k,k') = M I 



dh 



^q(k)^*(k') 



(27r)3 p2 



q^ -|- ie 



(7) 



when expressed in terms of the Coulomb wavefunctions for the center-of-mass energy 
E = 1? jM. As illustrated in Fig. 2 it includes the free propagator plus the effects of one. 
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Figure 2: The full Coulomb propagator is formed from the infinite sum of exchanged static 
photons. 



two and more exchanged static photons. In this way we find for the diagram in Fig. lb the 
value CqB{p) where 



B{p) = vsi^y 

is a convergent integral. 
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(27r)3 (2^)3 k2 + 72 



Gc{E-'k,\^) 



(8) 



The protons can suffer another rescattering as in Fig.lc before they are converted 
to a deuteron by the weak interaction. Denoting this extra bubble in the diagram by 
Jo{p), the full diagram is then just Co Jo times the previous diagram. The magnitude 
of the bubble equals the probability amplitude for the two protons to propagate from 
zero spatial separation back to zero separation, i.e. Jo(p) = Gc{E;0,0). In terms of the 
representation (0) of the Coulomb propagator it is given by the integral 



Mp) = M 



d^q 27rri{q) 



(2vr) 



3 p2-Kr)(q) 



Ip'^ — q^ + ie 



(9) 



which has an ultraviolet divergence. After regularization, it will give a renormalization of 
the coupling constant Cq. 

We can now continue to add in more such rescattering diagrams in Fig.lc. They form 
a geometric series which sums up to Co/(l — Co Jo)- All the diagrams including the first 
then gives for the hadronic part of the full transition amplitude 



Tf,{p) = A{p) + B{p) 



Co 



1 - Co Job) 



^p(O) 



(10) 



where the last factor tpp{0) = Crje^"^ gives the amplitude for the two incoming protons to 
meet at the first vertex. It can be written in a more compact and recognizable form by 
introducing the proton-proton scattering state 



\^p) = [l + Y.{GcVoT]\^, 



(11) 



n=l 



where the strong interaction potential (p|Vo|q) = Co is included to all orders in addition 
to the Coulomb interaction. Then we see that the matrix element (|^) is just the overlap 
integral between this wavefunction and the deuteron wavefunction (0), 
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d^k 

(2^ 



V'Kk)^p(k) 



(12) 



as follows from the expressions in (y) and ^ . This form of the transition matrix element 
was written down first by Bethe and Chritchfield[P and used subsequently by everyone 
considering the process in potential models. 



We can now evaluate the different parts of the transition matrix element (10). The first 
part (^ is most easily found in coordinate space where we have the Coulomb wavefunction 
(I). It gives 

/•OO 

A{p) = ^/^Cr,e''"' drre-^^+'P^''M{l-iri,2;2ipr) (13) 

Jo 



— . .0 2^1 i - «^, ^, ^, —T- 

(7 + ipY \ 7 + ipj 



Now the hypergeometric function 2Fi{0'-, b, b; z) = {1 — z) "■ so that the final result can be 
written as 



Aip) = C.e^'^o^^^e'^'^''^*'^"^?) (14) 

p + 7 

In the expression (P) for B{p) we notice that the integral over k' gives the complex con- 
jugate value of the Coulomb wavefunction at the origin. It therefore takes the form 

^ ' J (27r)3 J (27r)3 k^ + 72 p2 - q2 + ^g ^^^ ' 

The integral over k is just the previous result for A{q) so that 

^^' J (27r)3 g2 + ^2p2_q2^ ^g g2^„(9) _ 1 ^^^^ 

It is seen that when the momentum of the incoming protons is non-zero the integral yields 
a complex result. 

The infinite sum over proton bubble diagrams in Fig.l is just the proton-proton strong 
scattering amplitude modified by Coulomb corrections! 15]. It is given by the bubble inte- 
gral (^. In order to regularize it we use the special PDS scheme constructed by Kaplan, 
Savage and Wise for this effective theory ||lO|| . It is based on ordinary dimensional regular- 
ization around d = 3 dimensions. The difference lies in that poles in d = 2 dimensions are 
subtracted. This gives rise to terms which depend on the regularization point ^. In the 
present case we obtain with e = 3 — d 
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Here Ce = 0.5772 ... is Euler's constant and the function 



H{r]) = ^{iri) + ^r-- Hir]) (17) 

lirj 



The divergent 1/e piece will be absorbed in counterterms representing electromagnetic 
interactions at shorter scales. This replaces the bare coupling constant Cq with the renor- 
malized value Co(/x). Matching the calculated proton-proton scattering amplitude to the 
experimental one, we can determine this coupling constant in terms of the measured scat- 
tering length Up which gives the cross-section in the zero-energy limit 1 15], 
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1 f^V^ , 1 3 



The part of the scattering amplitude which is needed in (10) is therefore 



(18) 



1 M 



h aMH{r]) 



(19) 



This is in general complex because of the function (17). Its real part is the more phe- 
nomenologically relevant function h{r]) = He ip{irj) — In r] while the imaginary part is simply 
C^/2rj. At very low energies the real part h{rj) = 1/(12^7^) -|- 0(ry~^) dominates since the 
imaginary part is then exponentially small. These functions are well-known in the context 



of proton-proton scattering[18| 



For the very small proton momenta we have in the Sun, one evaluates the transition 
matrix element just as well at zero momentum Q. The first term (113) is then simply 



\A{p 



ISttC^ 



l^x 



ryO 



(20) 



where the parameter x = aM/^. It is therefore natural to introduce the standard reduced 
matrix element S 



A(p) 



ryO 



^■kC^ 
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Using 2'Krj[q) as a new integration variable in the expression (15) for B{p 
proportional to the integral 



i{x) 



°° , 2x 
dx- 



^^arctan(^) 



/o e-^ — 1 x^ + TT^x^ 
which we can only do numerically. We thus have the result 

A(0) = e^ -aMapIix) 



(21) 
0), it becomes 

(22) 
(23) 



The first part is identical with what one obtains in potential models Q while the second 
part seems to have a different dependence on the parameter x- However, by numerical 
integration, we find that it is in fact exactly the same. So far we have been unable to show 
this equality by analytical methods. 



However, it can be demonstrated by taking first the zero energy limit of the proton- 
proton wavefunction in the integrand of the transition matrix element ( |l2|) and then 
afterwards integrate. The regular Coulomb wavefunction (^) simplifies then to|16| 



a 



Foip) ~ -y^p'/'h{2y/2^) 



(24) 

when 2r] ^ p. The first part iMj of the matrix element will now be given by the integral 

/47r7 1 



A{p ^ 0) = C^e''^^ 



1] p' 



POO 

Jo 



which can be expressed in terms of a Whittaker function, 

|A(p^O)| = 



87rC2 x/2 

T X '2 



Now M_-^ i{x) = X^^ and thus we reproduce the first part of (23). 

In order to calculate B{p ^ 0) we go back to the result (^. Now we use a different 
representation of the Coulomb propagator constructed from the regular Fi^p) and the 
irregular Ge{p) eigenfunctionsO]. We only need the radial part in the 5-channel 



Gc{E;r,r') 



Mp Fo(p<)Go(p>) 
4vr p< p> 



(25) 



where r< and r> are the smallest and largest of r and r' respectively. The irregular 
solution Ge{p) is normalized so that the Wronskian GqFq — FqGq = 1. In (^) this Green's 
function is seen to enter with the argument r' = 0. Since Fq{p)/ p = Grj in the limit p — > 0, 
it then simplifies to 



Gc{E;r,0)e=o = -G 



MpGojp) 
Airp 



In the zero-energy limit 2i] ^ p we can then use 

Go{p) ^ ^VWpKi{2./2^p) 

For the function B{p) we thus find in this limit 

M [lyq 



(26) 



B{p -^ 0) 



/ dpp^/^e-^P/PKi{2^/2^) 
Jo 



P \ TT 

This integral is now given by the other Whittaker function, 



B{p -^ 0) 



An 



7^ X 



-w.Ux) 



Since W_-^ i (x) = W_^ _i(x) we can thus express the result by the simpler integral 

fOO ^ — t 
<X 



W^_i_i(x) = Xe^/'/ dt-^ (27) 

'2 Jy t 



By partial integration it can be written in terms of the exponential integral function Ei{x)- 
In this way we finally obtain for the original integral d2^ 



I{x) = --e''E^{x) (28) 

X 

in agreement with the standard result [^. 

With the previous value 7 = 45.71 MeV for the deuteron, we have x = 0-15 and 
the integral /(0.15) = 4.96. Combined with the measured value Up = — 7.82fm for the 
scattering length we then have A(0) = 2.51 for the reduced matrix element. A recent 
and most accurate calculation including higher order corrections based on a potential 
modelQ gives a value A^(0) = 7.05 ± 0.02. This corresponds to A(0) = 2.66. Our leading 
order result from effective field theory which we have found to agree with the zero-range 
approximation of potential models, is thus accurate to within 6% percent of the full result. 

In next order of the momentum expansion of the underlying effective field theory it is 
not clear if the transition matrix element can still be given as a simple overlap integral of 



the two wavefunctions as in (12). The corresponding Feynman diagrams for the transition 
matrix element are also technically more difficult to evaluate because they involve the 
Coulomb propagator in new ways. But it is important to calculate these corrections and 
verify that they are as small as is found in potential models. 
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